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We theoretically investigate the quantum interference of entangled two-photon states generated 
qq . in a nonlinear crystal pumped by femtosecond optical pulses. Interference patterns generated by the 

' polarization analog of the Hong-Ou-Mandel interferometer are studied. Attention is devoted to the 

^\ , effects of the pump-pulse profile (pulse duration and chirp) and the second-order dispersion in both 

the nonlinear crystal and the interferometer's optical elements. Dispersion causes the interference 
pattern to have an asymmetric shape. Dispersion cancellation occurs in some cases. 



o 

Q 

00 



I. INTRODUCTION 

Significant consideration has recently been given to the process of spontaneous parametric down-conversion in non- 



tj/y been 

linear crystals pumped by cw lasers |lf,||,|[Q] . The nonclassical properties of entangled two-photon light generated by 
this process have been used in many experimental schemes to elucidate distinctions between the predictions of classical 
^ ■ and quantum physics 0. Coincidence-count measurements with entangled two-photon states have revealed violations 
of Bell's inequalities and have been considered for use in nonclassical imaging and quantum cryptography Bj . 

A new frontier in these efforts is the generation of quantum states with three correlated particles (GHZ states) [p|,pO[, 
which would be most useful for further tests of the predictions of quantum mechanics. One way to create such states 
is to make use of pairs of two-photon entangled states that are synchronized in time, i.e., generated within a sharp 
time window EJ. This can be achieved by using femtosecond pump beams. Also, successful quantum teleportation 
! has already been observed using femtosecond pumping |l^] . 

For these reasons, the theoretical and experimental properties of pulsed spontaneous parametric down-conversion 
have been scrutinized |lj|,[^,[u|[l6| . It has been shown that ultrashort pumping leads to a loss of visibility of the 
coincidence-count interference pattern in type-II parametric down-conversion jh|,[l4|jl^] , and narrowband frequency 
03 filters are required to restore the visibility ]Tl| , |l^ , |l5| . 

This paper is devoted to a theoretical investigation of dispersion effects in femtosecond-pulsed spontaneous para- 
^ ' metric down-conversion. Particular attention is given to the effects of pump-pulse chirp and second-order dispersion 
, (in both the pump and down-converted beams) on the visibility and shape of the photon-coincidence pattern gener- 
ated by the polarization analog of the Hong-Ou-Mandel interferometer p7| . Dispersion cancellation, which has been 
extensively studied in the case of cw pumping ]l8| , is also predicted to occur under certain conditions for femtosecond 
■ down-converted pairs. 



II. SPONTANEOUS PARAMETRIC DOWN-CONVERSION WITH AN ULTRASHORT PUMP PULSE 

We consider a nonlinear crystal pumped by a strong coherent-state field. Nonlinear interaction then leads to the 
spontaneous generation of two down-converted fields (the signal and the idler) which are mutually strongly correlated 
0. Such a correlation can be conveniently described in terms of the two-photon amplitude A12 which is defined 
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as a matrix element of the product of electric-field operators E[ + \z\,ti) and {z 2l t 2 ) sandwiched between the 
entangled two-photon state \ip^) (for details, see Appendix A) and the vacuum state |vac): 

A 12 (z 1 ,t 1 ,z 2 ,t 2 ) = (v aC \E[ +) (z 1 ,t 1 )E 2 :+ \z 2 ,t 2 )\^(0,t)). (1) 
The positive-frequency part E^ of the electric-field operator of the jth beam is defined as 

E f \ z 3^j) = e j(%)/j( w fc 3 )%(%) exviikjZj - iuj kj t 3 ), j = 1, 2, (2) 

kj 

where stands for the annihilation operator of the mode with wave vector kj, ej(kj) denotes the normalization 
amplitude of the mode kj , and fj (wk ) characterizes an external frequency filter placed in the jth beam. The symbols 
k\ and fc 2 denote wave vectors in vacuum. 

At the termination of the nonlinear interaction in the crystal, the down-converted fields evolve according to free- 
field evolution and thus the two-photon amplitude A12 depends only on the differences t\ — t and t 2 — t. When the 
down-converted beams propagate through a dispersive material of the length I, the entangled two-photon state \ip^) 
given in Eq. (A4) in Appendix A provides the expression for Ai 2 j: 

,0 

-4i2,;( T i,T 2 ) = C / dz ^^/i(a) fel )^/ 2 (wfc 2 )fp +) (0,w fcp - ujp)cxp[i(k p - ki - k 2 )z] 

x cxp \i(k\ + fc 2 )^J 5(LUk p -tJk! - u k2 ) cxp[-«w fel Ti] CXp[-iWfc 2 T 2 ]. (3) 

The times T\ and t 2 are given as follows: 

-iu) kj Tj = ikpj - iuj kj {tj - t), j = 1 , 2. (4) 

The symbol £ p + \o, LUk p — uo^) denotes the positive-frequency part of the envelope of the pump-beam electric-field 
amplitude at the output plane of the crystal and cu® stands for the central frequency of the pump beam; the wave 

vectors k p , k\, and k 2 {k\ and k 2 ) are appropriate for the nonlinear crystal (dispersive material). The symbol L means 
the length of the crystal. The amplitudes ei(fci) and e 2 (fc 2 ) from Eq. (2) are absorbed into the constant C. 
A typical experimental setup for coincidence-count measurement is shown in Fig. 1. 
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FIG. 1. Sketch of the system under consideration: a pump pulse at the frequency lu p generates down-converted photons at 
frequencies coi and LO2 in the nonlinear crystal NLC. These waves propagate through a delay line of length I and are detected 
at the detectors Da and Db; BS denotes a beamsplitter; Aa and Ab are analyzers; Fa and Fb are frequency filters; and C 
indicates a coincidence device. 



We consider type-II parametric down-conversion for this exposition. In this case two mutually perpendicularly po- 
larized photons are provided at the output plane of the crystal. They propagate through a bircfringent material of a 
variable length I and then impinge on a 50/50 beamsplitter. Finally they are detected at the detectors and Db- 
The coincidence-count rate R c is measured by a coincidence device C. The beams might be filtered by the frequency 
filters F^ and F^ which can be placed in front of the detectors. Analyzers rotated by 45 degrees with respect to 
the ordinary and extraordinary axes of the nonlinear crystal enable quantum interference between two paths to be 
observed; either a photon from beam f is detected by the detector Da and a photon from beam 2 by the detector 
Db , or vice versa. 
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Including the effects of the beamsplitter and analyzers, the coincidence-count rate R c can be determined as follows 
0Q: 
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R*(l) = \ I 'II a I 'I'd \A V2 .,\l A .i u ) ■Avj.iHij.i S :! J ■ (•->) 
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where the two-photon amplitude A12.1 is given in Eq. (3). 

The normalized coincidence-count rate R n is then expressed in the form: 

R n {l) = l- P {l), (6) 

where 

P(1) = 7TB- dt A dt B Re[Ai2,i(tA,t B )A$ 2l (t B ,t A )], (?) 

ZK J -00 J-00 

and 

/>oo />oo 

R o=~ dt A dt B \A 12 ,i(tA,t B )\ 2 • (8) 

^ J — 00 J —oc 

The symbol Re denotes the real part of its argument. 

III. SPECIFIC MODELS INCLUDING SECOND-ORDER DISPERSION 

Let us assume that the nonlinear crystal and the optical material in the path of the down-converted photons are both 
dispersive. We proceed to generalize the models provided in Refs. |13],|l4|,[l5| by including the effects of second-order 
dispersion. 

The wave vectors k p (u)k p ), and /c2(wfc 2 ) of the beams in the nonlinear crystal can be expressed in the 

following form, when the effects of material dispersion up to the second order are included | fL9| : 

fcjKH^ + ^K - u O) + £L( Ukj _ u oy } i=p,l,2. (9) 



The inverse of group velocity 1/vj , and the second-order dispersion coefficient Dj , are given by 

1 dkj 



(10) 



Oi=2^U», i=P,l,2. (ID 



hi 



The symbol oj® denotes the central frequency of beam j. The wave vector k® is defined by the relation k® — kj(bj®). 

Similarly, the wave vectors k\ (co^ ) and k 2 {uJk 2 ) of the down-converted beams in a dispersive material outside the 
crystal can be expressed as: 



hi K ) = fc° + - K - w?) + ^ K - ^°) 2 ' i = 1 . 2 > ( 1 2 » 

9j ^ 

where 



1 dkj 



(13) 



and fc° = jfej(w?) 
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We further assume that frequency filters with a Gaussian profile, and centered around the central frequencies uj 
and oj 2 , are incorporated: 



fj(uk } ) = exp 



3 =1,2, 



(15) 



where cr* is the frequency width of the jth filter. 

Assuming frequency- and wave- vector phase matching for the central frequencies (uj p = uj® + ^2) an( l central wave 
vectors (k p = fcj + fcg), respectively, the two-photon amplitude *4i2,z(ti, t 2 ) defined in Eq. (3) can be expressed in the 
form: 



A 12 ,i(t 1 ,t 2 ) = C^exp(-iw5 l T 1 )exp(-i^r2) / efe / £ p +) (0, fl p ) 
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The frequencies flj, flj = uj]~ j — u)®, for j = 1, 2,p have been introduced in Eq. (16); Cj± denotes a constant. 

We proceed to devote further attention to special cases. We first consider an ultrashort pump pulse with a Gaussian 
profile: the envelope £p + \o, t) of the pump pulse at the output plane of the crystal then assumes the form p3]: 



£ p +\o,t) = Z p0 exp(- 



1 + ia 



-r 



(17) 



where £ p o is the amplitude, Td is the pulse duration, and the parameter a describes the chirp of the pulse. 
The complex spectrum £p + \z,fl p ) of the envelope £ p (z,t) is defined by 



For a pulse of the form given in Eq. (17) we obtain: 



1 f°° 

£(+\z,Sl p ) = — I dt£ p +\z,t)ex.p(iQ p t). 



4 +) (o,fi P ) = e f 
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(18) 



(19) 



where £ p = £ p o exp[— i arctan(a)/2] . 

Substituting Eq. (19) into Eq. (16) and using the identity 



/oc poo 
dQi / dfl 2 exp [— aiQ 2 — a 2 Q 2 — 20:12^1^2 + ioi^i — £1X2^2] 
-00 J — 00 
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we arrive at the following expression for the two-photon amplitude A\ 2 .i{ti 1 t 2 ): 



Al2,l(Ti,T 2 ) = C A 



■. exp(-iw5 ) ri) exp(-iw°T2)A 12i i(Ti, r 2 ) 



c 2 ^ 2 + c 2 2 Pi + 2 1C ic 2 



r 1 

^i2,i( T i,T 2 ) = / dz =exp 
J-l V P1P2 - 7 



4(/3i/3 2 - l 2 ) 



(20) 

(21) 
(22) 



The functions (3j(z), Cj(z), and 7(2) are defined as follows: 
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(3 j (z) = - s + b(l-ia)-i^l-i^ 
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7(2) =b(l-ia) -i^ z - 
The parameter & is a characteristic parameter of the pump pulse: 
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The quantities p(l) and i?o are then determined in accordance with their definitions in Eqs. (7) and (8), respectively. 
The quantity p(l) as a function of the length I of the bircfringent material then takes the form (wj = w° 1S assumed): 
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The functions /3j(2i,2 2 ), Cj(zi,2 2 ), and 7(21,22) are expressed as follows: 
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Similarly, the normalization constant i?o is given by the expression: 



2\c \2^2 r° 



Rq — 



2V1 + a 2 



C?29 



: exp 



7 2 



g 2 /3 2 + Clffl + gTClCg 

4(/3!/3 2 - 7 2 ) 



(26) 
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where 



^•(21, 22) = \ + 26 - Dj ( Zl - 22), j = 1, 2 

4^ 
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(28) 



It is convenient to consider the pump pulse characteristics at the output plane of the crystal, i.e., to use the 
parameters td and a. They can be expressed in terms of the parameters Toi and dj appropriate for the input plane 
of the crystal: 



a = 
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(29) 



In this case, the parameter bi 



'Di 



' 4(1 + a?) 

has the same value as the parameter 6 defined in Eq. (24). 
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Ignoring second-order dispersion in all modes [D p = D\ = D2 = 0), Eq. (25) reduces to the following analytical 
expression for the quantity p: 
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(31) 
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(32) 



Ar ; =ti- DL/2. 



(33) 



The symbol erf denotes the error function. When deriving Eq. (31) the condition D > was assumed. In Eq. (33), 
ti denotes the relative time delay of the down-converted beams in a birefringent material of length / and is defined as 
follows: 



n = , J- - J- M . 

,92 9i, 



(34) 



When second-order dispersion in the down-converted fields is omitted, the interference pattern can be determined 
for an arbitrary pump-pulse profile in terms of the autocorrelation function of the pump pulse. For details, see 
Appendix B. 



IV. DISCUSSION 



We now proceed to examine the behavior of the normalized coincidence-count rate R n on various parameters, from 
both analytical and numerical points of view. 

The profile of the interference dip in the coincidence-count rate Jllj (described by p as a function of I), formed by 
the overlap of a pair of two-photon amplitudes, can be understood as follows. The expression in Eq. (7) for p(l) can 
be rewritten in the form: 

^-KO J -00 J -00 



where 



t= A ~t tB , r = t A -t B , (36) 



and A r 12 1 — R- e [-^i2,z]; A\ 2 i = Im[„4 12i z]. The symbol Im denotes the imaginary part of the argument. Hence, 
according to Eq. (35), the overlaps of the real and imaginary parts of the two-photon amplitudes Ai2,i(t,r) and 
Ai2,i(t, —t) determine the values of the interference term p. The amplitude A\2,i{t, — r) can be considered as a mirror 
image of the amplitude -4i2,i(^ r ) with respect to the plane r = 0. When only first-order dispersion in the optical 
material is taken into account, the shape of the two-photon amplitude Ai2,i{t, r) does not depend on the length /; as 
I increases, the amplitude -4i2,z(i, r) moves only in the t-T plane. The shift in the r-direction is important, because 
it changes the degree of overlap of the amplitudes. This reveals the origin of the shape of the dip. 

The overlap of the two-photon amplitudes can be interpreted from the point-of-view of distinguishability of two 
paths leading to coincidence detection p3| . When the overlap is complete, the two paths cannot be distinguished 
and the interference pattern has maximum visibility. Incomplete overlap means that the paths can be "partially 
distinguished" and thus the visibility is reduced. 

We consider, in turn, the role played by pump-pulse duration and chirp, second-order dispersion in the nonlin- 
ear down-converting medium, second-order dispersion in the optical elements of the interferometer, and dispersion 
cancellation. 
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A. Pump-pulse duration and chirp 



In the absence of second-order dispersion and frequency filters, a useful analytical expression for the two-photon 
amplitude Ai 2: i=o(t, r) can be obtained: 

-—{ t + ° r )\- (37) 

The coefficients D and A are defined in Eq. (32). Equation (37) elucidates the role of pump-pulse parameters as 
discussed below. It is well known that for a cw-pump field the coincidence-count rate R n {ri) forms a triangular dip 
of width DL fly. The visibility is 100%, indicating maximum interference. An ultrashort pump pulse of duration roi 
leads to a loss of visibility (see Fig. 2) but the width of the dip remains unchanged This can be understood from 
the shape of the two-photon amplitude A\2^ = a{t, r) given in Eq. (37). In the r-direction the two-photon amplitude is 
confined to the region < r < DL for either cw or an ultrashort pump pulse; this confinement is responsible for the 
width of the dip. The two-photon amplitude is confined in the t-direction by the ultrashort pump-pulse duration [see 
Eq. (37)]. The tilt (given by the ratio A/D, sec Eq. (37)) of the amplitude in the t-r plane leads to a loss of visibility 
since the overlap of the amplitudes ^4i2,z(i, r) and Ai2j(t, — r) for a given optimum value of I cannot be complete for 
a nonzero tilt. The shorter the pump-pulse duration, the smaller the overlap, and the lower values of visibility that 
result. However, when values of the first-order dispersion parameters are chosen such that A = 0, the tilt is zero [see 
Eq. (37)] and no loss of visibility occurs as the pump-pulse duration shortens (for details, see pif), 
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FIG. 2. Visibility V (V = p/(2 — p)) as a function of the pump-pulse duration td%\ L — 3 mm, a = oo nm, and tij = 0; 
values of the other parameters are zero. In Figs. 2-8, the following parameters apply. Values of the inverse group velocities 
appropriate for the BBO crystal with type-II interaction at the pump wavelength X p — 397.5 nm, and at down-conversion 
wavelengths Ai = A 2 = 795 nm are: l/v p = 57.05 x KT 13 s/mm, 1/vi = 56.2 x 1CT 13 s/mm, and l/v 2 = 54.26 x KT 13 
s/mm. We assume that the optical materials for the interferometer are quartz, for which 1/gi = 51.81 x 10~ 13 s/mm and 
1/02 = 52.08 x 10~ 13 s/mm. 



Al2 ^ o(< ' r) = Tm \D\ rect Km) exp 




As indicated by the Eq. (37) for the amplitude ^12^=0, pump-pulse chirp (characterized by ai) introduces a phase 
modulation of the two-photon amplitude in the i-dircction. This modulation decreases the overall overlap of the 
corresponding two-photon amplitudes, given as a sum of the overlaps of their real and imaginary parts. Increasing 
values of the chirp parameter ai thus lead to a reduction of visibility. However, the width of the dip does not 
change. In fact, it is the parameter bi given in Eq. (30), combining both the pulse duration Tr>% and the chirp 
parameter a,, that determines the visibility in case of a Gaussian pump pulse. To be more specific the parameter 
bi is determined by the bandwidth Af2 p [AO p = V2\/l + af/r^i, see Eq. (19)] of the pump pulse according to the 
relation bi = l/[2(Ai7 p ) 2 ]. Thus, more generally, it is the bandwidth of the pump pulse that determines the interference 
pattern. As a consequence, dispersion of the pump beam between the pump-pulse source and the nonlinear crystal 
does not influence the interference pattern. 

Examination of Eqs. (B3) and (B4) in Appendix B shows that the dip remains symmetric since the function p(Ari) 
in Eq. (B3) is an odd function of At; for an arbitrary pump-pulse profile. 

Frequency filters inserted into the down-converted beams serve to broaden the two-photon amplitude Ai2,;(t, t) 
both in the t- and r-direction. Broadening in the r-direction leads to wider dips, whereas that in the t-direction 
smooths out the effect of tilt discussed above and thereby results in a higher visibility. The narrower the spectrum 
of frequency filters, the wider the dip, and the higher the observed visibility. The effect of chirp is suppressed by the 
presence of frequency filters, because they effectively make the complex pump-pulse spectrum narrower and hence 
diminish relative phase changes across such a narrowed complex spectrum. 
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B. Second-order dispersion in the nonlinear crystal 



Second-order dispersion in the pump beam causes changes in the pulse phase (chirp) as the pulse propagates and 
this leads to broadening of the pulse. The effect of such pump-pulse broadening is transferred to the down-converted 
beams, as is clearly shown by the behavior of the two-photon amplitude Ai2.i(t,r) illustrated in Fig. 3. 




FIG. 3. Absolute value of the two-photon amplitude A — \Ai2,i=o{t, r)| for nonzero second-order dispersion of the pump 
beam; the variables t and r are in units of 1CF 13 s; td% = 1-55 x ICP 13 s, L = 3 mm, a = 100 nm, D p = 1 x 10~ 25 s 2 /mm, and 
di = 0; values of the other parameters are zero. 

In this figure, the amplitude in the region near t = 0s has its origin near the output plane of the crystal where the 
pump pulse is already broadened as a result of its having propagated through the dispersive crystal. At the other 
edge, near t w 6 x 10~ 13 s the down-converted light arises from the beginning of the crystal where the pump pulse 
has not yet suffered dispersive broadening. The profile of the interference dip is modified as follows: An increase in 
the second-order dispersion parameter D p leads to an increase of visibility, but no change in the width of the dip, as 
illustrated in Fig. 4(a). For appropriately chosen values of D p a small local peak emerges at the bottom of the dip [see 
Fig. 4(a)]. Nonzero initial chirp (a^) of the pump beam can provide a higher central peak but on the other hand it 
reduces the visibility [see Fig. 4(b)]. The peak remains, but is suppressed, in the presence of narrow frequency filters. 




0.4 0.6 

-5-3-1135 -5-3-1135 



(a) At, (1CT 13 s) ( b ) At, (1CT 13 s) 

FIG. 4. Coincidence-count rate -R„(At;) (a) for various values of the second-order dispersion parameter D p : D p = s 2 /mm 
(plain curve), D p = 5 x 10~ 26 s 2 /mm (*), D p = 1 x 10~ 25 s 2 /mm (A), and C„ = 3x 10~ 25 s 2 /mm (O), en = and (b) for 
various values of the chirp parameter en: en — (plain curve) and <n = 2 (*), D p — 5 x 10 -26 s 2 /mm; td% = 1.55 x 10~ 13 s; 
L = 3 mm; a — 50 nm; values of the other parameters are zero. 

Now we turn to second-order dispersion in the down- converted beams (nonzero D\, D2), which broadens the two- 
photon amplitude A 12j ;(i, r) in the r-, as well as in the i-direction. As demonstrated in Fig. 5, this leads to a 
broadening of the dip, as well as asymmetry and oscillations at its borders. When values of D\ increase, visibility 
decreases at first and then later increases. Nonzero chirp leads to a lower visibility, but tends to suppress ocillations 
at the borders of the dip. Frequency filters, which behave as discussed above, suppress asymmetry 
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FIG. 5. Coincidence-count rate R n (Ari) in case of second-order dispersion in beam 1 (plain curve, D\ = 1 x 1CP 25 s 2 /mm, 
D 2 = s 2 /mm) and in beam 2 (*, Di = s 2 /mm, D 2 = 1 x 10~ 25 s 2 /mm); r Di = 1.55 x 10~ 13 s, L — 3 mm, and <7i = cr 2 = 50 
nm; values of the other parameters are zero. 

When second-order dispersion occurs in all three modes, the two-photon amplitude A\2,i(t, r) is broadened for 
smaller values of r (mainly owing to dispersion in the pump beam) as well as for greater values of r (mainly owing to 
dispersion in the down-converted beams). As a result, the interference pattern comprises all of the features discussed 
above: a local peak may emerge at the bottom of the dip, the dip is broadened and asymmetric, and there occur 
oscillations at the borders of the dip. 

To observe the above mentioned effects caused by dispersion in a nonlinear crystal, relatively large values of the 
dispersion parameters D p , Di, and D 2 are required. For example, our simulations make use of parameter values 
that are approximately an order of magnitude higher than those of the BBO crystals commonly used in type-II 
down-conversion-based interferometric experiments. 

C. Second-order dispersion in the interferometer's optical elements 

Second-order dispersion in an optical material (di, d 2 ) through which down-converted photons propagate leads to 
asymmetry of the dip. The dip is particularly stretched to larger values of I (see Fig. 6) as a consequence of the 
deformation and lengthtcning of the two-photon amplitude Ayi.i in a dispersive material. The higher the difference 
d\ —d 2 of the dispersion parameters, the higher the asymmetry and the wider the dip; moreover its minimum is shifted 
further to smaller values of I (see Fig. 6). Asymmetry of the dip is also preserved when relatively narrow frequency 
filters are used though the narrowest filters remove it. Chirp decreases visibility but the shape of the dip remains 
unchanged. 




0.5 -I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

10 20 30 

1 (mm) 

FIG. 6. Coincidence-count rate R n (l) for various values of the second-order dispersion parameter d — di — d 2 of an optical 
material; d = s 2 /mm (plain curve), d = 1 x 1CT 26 s 2 /mm (*), d — 5 x 1CT 26 s 2 /mm (A), and d = 1 x 1CT 25 s 2 /mm (O); 
tdi = 1-55 x 1CF 13 s, L — 3 mm, and o\ = <T2 = 50 nm; values of the other parameters are zero. 



D. Dispersion cancellation 

Asymmetry of the dip caused by second-order dispersion in an optical material through which down-converted 
photons propagate can be suppressed in two cases. In the first case, for a pump pulse of arbitrary duration, dispersion 
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cancellation occurs when the magnitude of second-order dispersion in the path of the first photon (given by d\l) equals 
that of the second photon (given by d^l). This observation immediately follows from Eqs. (25) and (26), in which 
the effect of second-order dispersion is prescribed by the parameter (d\ — <i 2 )/. Dispersion cancellation is a result 
of completely destructive interference between the amplitudes A\2,i{t,r) and A\2,i{t, — t) for which there is nonzero 
overlap. This is demonstrated in Fig. 7 for I = 25 mm, i.e. for which p = 0. 




FIG. 7. Absolute value of the two-photon amplitude A — |Ai2,i(i, t)| for the same amount of second-order dispersion in 
the down-converted beams (di = di = 1 x 1CP 25 s 2 /mm) for I = 25 mm; the variables t and r are in units of 1CP 13 s; 
Tm = 1.55 x 1CT 13 s, L — 3 mm, and a = 100 nm; values of the other parameters are zero. 

When the pulse duration is sufficiently long (in the cw regime) dispersion cancellation occurs for arbitrary mag- 
nitudes of second-order dispersion (given by d\l and efeO present in the paths of the down-converted photons. The 
gradual suppression of the asymmetry of the dip as the pump pulse duration increases is shown in Fig. 8. 




0.0 A ' 1 1 r 1 1 1 1 1 1 1 1 

-10 10 20 30 40 50 

1 (mm) 



FIG. 8. Coincidence-count rate R„(l) shows a gradual suppression of dispersion effects (asymmetry) as the pump-pulse dura- 
tion increases; r Di = 1.55 x 10~ 13 s (plain solid curve), r Di = 5x 10~ 13 s (*), and r Di = 1 x 10~ n s (O), d — di —di = 5 x 10~ 26 
s 2 /mm; for comparison tch = 1-55 x 10~ 13 s, d = s 2 /mm (dashed curve); L = 3 mm; o\ = = 50 nm; values of the other 
parameters are zero. 

Dispersion cancellation has its origin in the entanglement of the photons, i.e., in the fact that the permitted values 
of the frequency ui\ and the frequency w 2 are governed by the relation 5(u) p — lo\ — uj 2 ) 1 where ui p lies within the 
pump-pulse spectrum. 

V. CONCLUSION 

We have developed a description of two-photon type-II spontaneous parametric down-conversion produced when 
ultrashort pulses from a femtosecond laser are used to pump an appropriate nonlinear medium, as well as the associated 
two-photon interference effects. The model includes frequency modulation of the pump pulse (chirp) and dispersion 
in both the nonlinear crystal and the interferometer's optical elements. The influence of these features on the depth 
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and asymmetry characteristics of a photon-coincidence interference dip have been established. 

We showed that the interference pattern is determined by the bandwidth of the pump pulse; the larger the band- 
width, the lower the interference-pattern visibility. This implies that dispersion of the pump beam before the nonlinear 
crystal does not influence the interference pattern. Second-order dispersion of the pump beam in the nonlinear crystal 
can result in the occurrence of a local peak at the bottom of the interference dip. Second-order dispersion of the 
down-converted photons in the crystal can result in oscillations at the borders of the dip, whereas dispersion of the 
down-covcrtcd photons in the interferometer's optical materials (e.g., the delay line) can produce an asymmetry in 
the dip. These effects can be used to measure the dispersion parameters of both a nonlinear crystal and an arbitrary 
optical material. Dispersion cancellation has been revealed for pump pulses of arbitrary duration when the amount 
of dispersion in the two down-converted beams is identical and in general for sufficiently long pump pulses. 
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APPENDIX A: DETERMINATION OF AN ENTANGLED TWO-PHOTON STATE 

The interaction Hamiltonian of the process of spontaneous parametric down-conversion can be written in the form 
Q: 

Hm(t)= / dz X {2) E( + \z,t)E[-\z,t)E { -\z,t)+h.c., (Al) 



—L 



where x*- 2 ^ is the second-order susceptibility, E p denotes the positive-frequency part of the electric-field amplitude 

of the pump field, and e[ ^ (E 2 ) is the negative-frequency part of the electric- field operator of down-converted field 
1 (2). The nonlinear crystal extends from z = —L to z = 0. The symbol h.c. means Hermitian conjugate. 

Expanding the interacting fields into harmonic plane waves, the interaction Hamiltonian iJj nt in Eq. (Al) can be 
recast into the form: 



H int (t) = C int [ dzJ2J2J2^ 2) 4 + HQ^k p -uj p )4(h)aUk2) 



^ k p k\ k-2 

x exp [i(k p - ki - k 2 )z - i(uj kp - u kl - u ks )t] + h.c, (A2) 



where C ln t is a constant. The symbol £ P +S> (0, uj kp — Lo p ) denotes the positive-frequency part of the envelope of the 
pump-beam electric-field amplitude at the output plane of the crystal; k p stands for the wave vector of a mode in the 
pump beam, and uj p stands for the central frequency of the pump beam. The symbol a\(ki) (62(^2)) represents the 
creation operator of the mode with wave vector k\ (£2) and frequency oj kl (w/c 2 ) m the down-converted field 1 (2). We 
note that the phases of all three interacting fields in space are chosen in such a way that they are zero at the output 
plane of the crystal. 

The wave function \ip( 2 \0,t)) describing an entangled two-photon state whose phases are set equal to at z = is 
given by: 



|^ (2) (0,i)) =4 I dt'H- mt (t')\™c), (A3) 



where |vac) denotes a multimode vacuum state. 

For times t sufficiently long so that the nonlinear interaction is complete, the entangled two-photon state \ip^ 2 \0, t)) 
can be obtained in the form: 

|</> (2) (0, t)) =C^J dzJ2J2J2 4 +) (°. - ( fc i)4(fc 2 ) exp [i(k p -kt- k 2 )z] 

k p k\ k 2 

x S{uj kp - Wfci - uj k2 ) exp [i(oJ kl + Uk 2 )t] |vac). (A4) 
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The susceptibility \^ 1S included in the constant C$. We note that for times during which the down-converted fields 
are being created in the crystal, the appropriate wave function differs from that in Eq. (A4). However, detectors 
are placed at a sufficiently large distance from the output plane of the crystal to assure that such "partially evolved" 
states cannot be detected. 



APPENDIX B: INTERFERENCE PATTERN FOR AN ARBITRARY PUMP-PULSE PROFILE 

We assume an arbitrary complex spectrum £p + \—L, O p ) for the envelope of the pump pulse at the input plane of 
the crystal. We further take into account the effect of second-order dispersion only in the pump beam and assume 
frequency filters of the same width (cti =02). Under these conditions, the normalized coincidence-count rate R n in 
Eq. (6) can be expressed in terms of the autocorrelation function of the pump field. 

?(+), 



Let us introduce the field £ pa '(z, t) according to the definition: 

D p (z + L) „ 2 
1 47T " p 



/OC 
-00 



cxp 



n 2 

p 



cxp(-iflpt), 



(Bl) 



where a = \/2o\. The above expression describes the propagation of the pump beam through a dispersive material 
(a multiplicative term describing first-order dispersion is not explicitly included here). Equation (Bl) also includes 
frequency filtering having its origin in the filtering of the down-converted beams and their entanglement with the 
pump beam. 

The two-photon amplitude Ai2, n (ti, T2) can then be derived from the expression in Eq. (16): 



Ai2, Ti (ti,t 2 ) = ^exp(-iwiVi)exp(-iw2T 2 ) 
,0 

V^t / dz£ { +\z, (n + ti + r 2 )/2 - Az) exp 



a 2 
-^(ti+t ; -t 2 + Dz) 



(B2) 



where the parameters D and A are defined in Eq. (32) and the relative time delay 77 of the down-converted beams is 
introduced in Eq. (34). 

The quantity p given in Eq. (7) then has the form (again it is assumed that uo\ — w"): 



p(Atz) 



|CU|V2mr<7 



x Re 



4i?o 

L/2 
-L/2 



dz\ 



L/2 



L/2 



dz2'ja(zi,Z2,A(zi - Z2)) exp 



An + —{zi + z 2 ) 



(B3) 



where Ar ; is defined in Eq. (33). The correlation function j a (zi,z 2 ,x) of two pulsed fields at positions Z\ and z 2 is 
written as 



/oo 
dt 4+> (zi - L/2, t)£t~\z2 - L/2, t + x). 
-OO 



(B4) 



The constant Ro occurring in Eq. (B3) is expressed as follows: 



Ro = 



|CU| 2 \/2W f L ' 2 



/L/2 /-L/2 
dzi / dz2ia(zi,Z2, A(zi - z 2 )) exp 
-L/2 J-L/2 



32 



-(z 1 -z 2 f 



(B5) 



For a Gaussian pulse with the complex spectrum as given in Eq. (19), the correlation function j a becomes 



Ja(zi, Z 2 ,A(z 1 - Z 2 )) 



\tp\ 



2V1 + o? VW^2) 



exp 



A 2 (zi - Z2) : 
Ail>(zi,z 2 ) 



2 D 

ip(Z!,Z 2 ) = 2bi + -j - i-r^(z! - z 2 ) 



(B6) 
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which, together with Eqs. (B3) and (B5), leads to expressions which agree with those derived from Eqs. (25) and 
(27). The parameter bi is defined in Eq. (30). 

The experimental setup without frequency filters (<r — > oo) is of particular interest. In this case, using the identity 
Y / 7rcrexp(— cr 2 y 2 /4) — > 2i:8(y) for a — > oo, Eqs. (B3) and (B5) provide a useful expression for the function p(An): 

7oo(U,0,(J)^ |y-L/2 

x 7oo (z,-z-2At ; /^,2z + 2At / /D)}, (B7) 
where rect(x) is the rectangular function (rect(x) = 1 for < x < 1 and rect(x) = otherwise). 
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